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Abstract 

Batteries are used in a wide variety of applications. In recent 
years, they have become popular as a source of power for 
electric vehicles such as cars, unmanned aerial vehicles, and 
commericial passenger aircraft. In such application domains, 
it becomes crucial to both monitor battery health and perfor- 
mance and to predict end of discharge (EOD) and end of use- 
ful life (EOL) events. To implement such technologies, it is 
crucial to understand how batteries work and to capture that 
knowledge in the form of models that can be used by moni- 
toring, diagnosis, and prognosis algorithms. In this work, we 
develop electrochemistry-based models of lithium-ion batter- 
ies that capture the significant electrochemical processes, are 
computationally efficient, capture the effects of aging, and are 
of suitable accuracy for reliable EOD prediction in a variety 
of usage profiles. This paper reports on the progress of such 
a model, with results demonstrating the model validity and 
accurate EOD predictions. 

1. Introduction 

With electric unmanned aerial vehicles, electric/hybrid cars, 
and commercial passenger aircraft, we are witnessing a dras- 
tic increase in the usage of batteries to power vehicles. These 
vehicles are being deployed in military, civilian and scien- 
tific applications all over the world. However, for battery- 
powered vehicles to operate at maximum efficiency and reli- 
ability, it becomes crucial to both monitor battery health and 
performance and to predict end of discharge (EOD) and end 
of useful life (EOL) events. 

For example, NASA’s Mars Global Surveyor stopped oper- 
ating in late 2006 due to battery overheating caused by the 
raditor being oriented directly towards the sun, resulting in a 
significant loss of capacity (Saha & Goebel, 2009). Both the 
AFRL ARGOS satellite and the Viking 2 Mars Lander also 
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suffered battery failures. More recently, problems have arisen 
in the Boeing Dreamliner aircraft, which utilize lithium-ion 
(Li-ion) batteries for auxiliary power (Ross, 2013). Bat- 
tery monitoring failed to capture thermal runaway conditions, 
leading to fire. With advanced health monitoring and predic- 
tive capabilities, such problems can be prevented. 

To implement such technologies, it is crucial to understand 
how batteries work and to capture that knowledge in the 
form of models that can be used by monitoring, diagnosis, 
and prognosis algorithms. Several battery modeling method- 
ologies have been developed in the literature. The models 
can generally be classified as empirical models, electrochem- 
ical engineering models, multi-physics models, and molec- 
ular/atomist models (Ramadesigan et al., 2012). Emprical 
models (e.g., (Saha & Goebel, 2009)), are based on fitting cer- 
tain functions to past experimental data, without making use 
of any physicochemical principles. Electrical circuit equiv- 
alent models are popular forms of empirical models, e.g., 
(Ceraolo, 2000; Chen & Rincon-Mora, 2006; Daigle, Sax- 
ena, & Goebel, 2012). Electrochemical engineering models 
are typically continuum models that include electrochemical 
kinetics and transport phenomena, e.g. (Karthikeyan, Sikha, 
& White, 2008; Rong & Pedram, 2006; Ning & Popov, 2004; 
Doyle, Fuller, & Newman, 1993; Newman & Tiedemann, 
1975). The former type of model has the advantage of be- 
ing computationally efficient, but the disadvantage of lim- 
ited accuracy and robustness, due to the approximations in 
battery behavior that are made, and as a result of such ap- 
proximations, cannot represent aging well. The latter type of 
model has the advantage of being very accurate, but are often 
computationally inefficient, having to solve complex sets of 
partial differential equations, and thus being ill-suited to on- 
line application of prognostics. Similarly, multi-physics (Lee, 
Smith, Pesaran, & Kim, 2013) and atomist models are even 
less suited to online application. 

In this work, we develop an electrochemistry-based model of 
Li-ion batteries, that (i) capture the significant electrochemi- 
cal processes, (ii) are computationally efficient, (Hi) capture 
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the effects of aging, and (iv) are of suitable accuracy for reli- 
able EOD prediction in a variety of usage profiles. The model 
can be considered an electrochemical engineering model, but 
unlike most such models found in the literature, we make cer- 
tain approximations that allow us to retain computational effi- 
ciency. Although we focus here on Li-ion batteries, the model 
is quite general and can be applied to different chemistries 
through a change of model parameter values. This paper re- 
ports on the progress of such a model, providing model vali- 
dation results and EOD prediction results. 

The paper is organized as follows. Section 2 provides 
background material on batteries. Section 3 describes the 
model-based prognostics framework. Section 4 develops the 
electrochemistry-based battery model. Section 5 provides 
model validation and battery prognosis results. Section 6 con- 
cludes the paper. 


2. Background 

A battery converts chemical energy into electrical energy, and 
often consists of many cells. A cell consists of a positive elec- 
trode and a negative electrode with electrolyte in which the 
ions can migrate. For Li-ion, a common chemistry is a pos- 
itive electrode consisting of lithium cobalt oxide (LLC 0 O 2 ) 
and negative electrode of lithiated carbon (Li E C). These ac- 
tive materials are bonded to metal-foil current collectors at 
both ends of the cell and electrically isolated by a microp- 
orous polymer separator film that is permeable to Li ions. The 
electrolyte enables lithium ions (Li + ) to diffuse between the 
positive and negative electrodes. The lithium ions insert or 
deinsert from the active material (known as intercalation and 
deintercalation) depending upon the electrode and whether 
the active process is charging or discharging. 

A schematic of the battery during the discharge process is 
shown in Fig. 1. When fully charged, the active (mobile) Li 
ions reside in the negative electrode. Upon connecting a load 
to the battery, current is allowed to flow from the positive 
to the negative electrode. This is supported by the oxidation 
reaction (loss of electrons) in the negative electrode (acting 
as the anode): 


Li x C i^g !>c + ;cLi++a;e - ) 


‘(t) Load 



+ Current Collector - Current Collector 

Figure 1. Li-ion battery during discharge. 


i(t) Current Source 



+ Current Collector - Current Collector 

Figure 2. Li-ion battery during charge. 


the negative to the positive electrode, as shown in Fig. 2. The 
active material in the positive electrode (acting as the anode) 
is oxidized and Li ions are de-intercalated with the corre- 
sponding reaction: 

LiCo0 2 > Li, _ ^.CoOa + xLi + + xe~ 

that results in the loss of Li ions and electrons, which can then 
move to the negative electrode (acting as the cathode). The 
reduction reaction takes place in the negative electrode: 


which results in the liberation of Li ions and electrons. The 
electrons move through the load from the negative to positive 
electrode, and the Li ions move in the same direction from 
the negative to the positive electrode through the separator. 
In the positive electrode (acting as the cathode), the reduction 
reaction (gain of electrons) takes place: 

Li, _ x Co0 2 + xLi+ + xe~ - ischarg - e > LiCo0 2 
During charging, a current source forces current to move from 


C + xLi + + xe~ Li^C 

Theoretically the cell has a voltage of around 4.2 V when 
fully charged. The terminal voltage of the battery rises/falls 
with a charge/discharge cycle, respectively. At the end of 
the charging/discharging period the battery voltage settles to 
a steady-state value that is a function of its state of charge 
(SOC). SOC is conventionally defined to be 1 when the bat- 
tery is fully charged and 0 when the battery is fully dis- 
charged. Hence even in a degraded battery the SOC for max- 
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imum charge is always 1 by definition. 

The charge/discharge dynamics of batteries can be character- 
ized by taking voltage measurements under constant current 
charge/discharge conditions. The rate of charge/discharge is 
measured relative to the battery capacity C. For example, a 

0. 1. discharge rate for a 5 Ah battery is 0.5 A. 

Charging typically consists of periods of constant current 
(CC) and/or constant voltage (CV) charging. At lower SOC, 
a CC charge is applied to bring the battery voltage to the CV 
level. Once the desired voltage level is achieved, the charger 
switches to CV mode and current decreases until 100% SOC 
is reached. The charger is not operated in CV mode at lower 
SOC levels since this will increase the amount of current 
flowing throught the battery while charging, leading to ad- 
ditional heat generation and side reactions taking place that 
may affect the battery life. The charge acceptance is the max- 
imum permissible rate at which batteries can be charged. For 
Li-ion batteries the recommended charge acceptance is C/3. 

The end of life of a battery depends on the chemistry, 
discharge-charge cycling, temperature, and storage condi- 
tions, among other factors. Atmospheric temperature affects 
the operating performance of the batteries. At low temper- 
atures, ionic diffusion and migration could be hindered and 
also damaging side reactions like lithium plating may take 
place. High temperatures cause corrosion and generation of 
gases leading to an increase in internal pressure. As the bat- 
tery ages, degradation results in the decrease in the observed 
battery capacity. This is primarily due to a loss of mobile ions 
due to parasitic or side reactions and an increase in the inter- 
nal resistance. Internal resistance leads to ohmic losses that 
generate heat and accelerate the aging process. Some relevant 
physical aging mechanisms in the electrodes are: 

1. Solid-electrolyte interface (SEI) layer growth: The neg- 
ative electrode degrades with the growth of the SEI layer 
leading to an increase in the impedance. The layers are 
formed during cycling and storage at high temperatures 
and entrains the lithium. 

2. Lithium corrosion: Lithium in the active carbon material 
of the negative electrode corrodes over time leading to 
degradation. This causes a decrease in the capacity due 
to irreversible loss of mobile lithium ions. 

3. Lithium plating: At low temperatures, high charge rates 
and low cell voltages forms a plating layer on the nega- 
tive electrode which leads to irreversible loss of lithium. 

4. Contact loss: SEI layer disconnects from the negative 
electrode which leads to contact loss and an increase in 
impedance. 

3. Battery Prognostics 

In this section we describe the general model-based prog- 
nostics problem and an implementation architecture, and de- 


scribe how it applies to batteries. 

3.1. Problem Formulation 

We assume the system model may be generally defined as 

x(fc + 1) = f (k,x.(k), 9(k),u(k),v(k)), (1) 

y{k) = h (k, x(fc), d(k), u (k), n (&)), (2) 

where k is the discrete time variable, x(fc) £ R”* is the 
state vector, 9(k) £ R™ 9 is the unknown parameter vector, 
u (k) £ R"“ is the input vector, v(fc) £ R"” is the process 
noise vector, f is the state equation, y(fc) £ R"« is the output 
vector, n (k) £ R"” 1 2 3 4 is the measurement noise vector, and h 
is the output equation. 1 

In prognostics, we are interested in predicting the occurrence 
of some event E that is defined with respect to the states, 
parameters, and inputs of the system. We define the event 
as the earliest instant that some event threshold Te : R" x x 
R" 9 x R n “ -A B, where B = {0, 1} changes from the value 
0 to 1. That is, the time of the event kE at some time of 
prediction kp is defined as 

k E (k P ) = inf{fc G N: k > k P A T E (x{k), 9{k), u(k)) = 1}. 

(3) 

The time remaining until that event, A k E , is defined as 

Ak E (kp) A k E (kp) — kp. (4) 

In the context of systems health management, T E is defined 
via a set of performance constraints that define what the ac- 
ceptable states of the system are, based on x(fc), 9(k), and 
u (k) (Daigle & Goebel, 2013). For batteries, there are two 
types of predictions that we are interested in. The first is 
EOD, defined by a voltage threshold V E od- In this case, T E 
is defined by V < V E o /)• that is, when the battery voltage 
is less than the cutoff voltage, EOD is reached and T E eval- 
uates to 1. The second type of prediction is EOL, which is 
typically defined by a lower bound on the effective battery 
capacity (Saha & Goebel, 2009). In this case, T E is defined 
by C < C E ol , where C is the measured battery capacity and 
Ceol is the lower bound on capacity. 

3.2. Prognostics Architecture 

We adopt a model-based prognostics architecture (Daigle 
& Goebel, 2013), in which there are two sequential prob- 
lems, (0 the estimation problem, which requires determin- 
ing a joint state-parameter estimate p(x(fc), 9(k)\y(ko'.k)) 
based on the history of observations up to time k, y{ko-k), 
and (ii) the prediction problem, which determines at 
kp, using p(x(k),9(k)\y(k(,:k)), a probability distribution 


'Here, we use bold typeface to denote vectors, and use n a to denote the 
length of a vector a. 
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Prognostics 



Figure 3. Prognostics architecture. 


p(kE(kp)\y(ko:kp)). The distribution for A fcg can be triv- 
ially computed from p(kE(kp)\y(ko'-kp)) by subtracting kp. 
Both steps require a dynamic model in the form of equations 1 
and 2. 

The prognostics architecture is shown in Fig. 3. In discrete 
time k, the system is provided with inputs it/, and provides 
measured outputs y/.. The estimation module uses this infor- 
mation, along with the system model, to compute an estimate 
p(x(fc), 0{k)\y(ko'.k)). The prediction module uses the joint 
state -parameter distribution and the system model, along with 
hypothesized future inputs, to compute the probability distri- 
bution p( k p(kp')\y( A'o :/>: /> ) ) at given prediction times kp. 

In the application of prognostics in this example, we do not 
account for any uncertainty except for that provided in the 
state estimate, as our focus is on determining how accurate 
the model can predict EOD given precise information about 
the future. A more general prognostics architecure that ac- 
counts for these additional sources of uncertainty is described 
in (Daigle, Saxena, & Goebel, 2012; Sankararaman, Daigle, 
Saxena, & Goebel, 2013). 

4. Battery Modeling 

In order to predict end-of-discharge as defined by a voltage 
cutoff, the battery model must compute the voltage as a func- 
tion of time given the current drawn from the battery. There 
are several electrochemical processes that contribute to the 
cell’s potential that make this a difficult problem. We fo- 
cus here on a lumped-parameter ordinary differential equa- 
tions form, so it is efficient and usable for on-line prognos- 
tics, yet still considers the main electrochemical processes. 
We focus here specifically on Li-ion 18650 batteries with an 
average nominal voltage of 3.7V and nominal capacity of 
2200mAh, however, the model is still general enough that 
with some modifications it may be applied to different bat- 
tery chemistries. 

The voltages of a battery are summarized in Fig. 4 (adapted 
from (Rahn & Wang, 2013)). The overall battery voltage 
V(t) is the difference between the potential at the positive 
current collector, (f> s ( 0, t), and the negative current collector, 
<t> s (L, t), minus resistance losses at the current collectors (not 
shown in the diagram). As shown in the figure, the potentials 
vary with the distance d, £ [0, L\, because the loss varies with 
distance from the current collectors. 





Figure 4. Battery voltages. 


The potentials at the current collectors are described by sev- 
eral voltage terms. At the positive current collector is the 
equilibrium potential Vjj , p . This voltage is then reduced by 
V SJ) , due to the solid-phase ohmic resistance, and V VtP , the 
surface overpotential. The electrolyte ohmic resistance then 
causes another drop V e . At the negative electrode, there is a 
drop V / ', / rj due to the surface overpotential, and a drop V s , n 
due to the solid-phase resistance. The voltage drops again 
due to the equilibrium potential at the negative current col- 
lector V[j n . We describe each of these voltage terms in turn. 


4.1. Equilibrium Potential 

In an ideal battery, the equilibrium potential of an electrode 
is described by the Nernst equation: 


RT 

Vu,i — Uo H — In 

nb 


'y fi ,i% [3 ,i 


(5) 


or, equivalently. 


RT 





7 / 3 , A 
7a, i) 


(6) 


where i refers to the electrode (n for negative or p for posi- 
tive), Uq is a reference potential, R is the universal gas con- 
stant, T is the electrode temperature, n is the number of elec- 
trons transferred in the reaction (n = 1 for Li-ion), F is Fara- 
day’s constant, x is the mole fraction, where a refers to the 
lithium-intercalated host material and /3 to the unoccupied 
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host material, and 7 is the activity coefficient (Karthikeyan 
et al., 2008). When discharging, Li ions move out of the 
negative electrode and into the positive electrode, so x a , n 
and x p p decrease while xp, n and x ap increase. As a result, 
Vjj, P — Vjj, n will decrease. In an ideal condition the activity 
coefficients are unity and so the last term disappears, resulting 
in the classical form of the Nernst equation. 

In the case of a Li-ion battery, xp p = 1 — x a> i, so we can 
rewrite the equation making this substitution. We then define 
x, as x a q and drop the a subscript: 

Vu,i = Uq H — In f — ^ + Vjf^r i 5 (7) 

nb \ Xi ) 

where Vint,* is the activity correction term (0 in the ideal con- 
dition). For Li-cCoC^, x must be at least 0.4; Li cannot be 
reversibly removed beyond that (Karthikeyan et ah, 2008). 

We let qi represent the amount of Li ions in electrode i, as 
measured in Coulombs. Recall that the flow of Li ions moves 
opposite to the flow of current. Therefore q, changes in the 
same direction as 27 , and we may define x, based on q, : using 


where < 7 max = q p + q n refers to the total amount of avail- 
able (i.e., mobile) Li ions. It follows then that x p + x n = 1. 
When fully charged, x p = 0.4 and x n = 0.6. When fully 
discharged, x p = 1 and x n = 0 . 

Batteries are nonideal and therefore the assumption of unity 
activity coefficients is not applicable. As a result, battery 
models often resort to curve-fitting for modeling the equilib- 
rium potential. The equilibrium potential can be obtained by 
discharging the battery at a very slow rate, such that a con- 
centration gradient will not build up, temperature is steady, 
and voltage drops due to internal resistances and other over- 
potentials are negligible. Fig. 5 shows the measured voltage 
(Vu p — V U n ) as a function of mole fraction against that pre- 
dicted by the Nernst equation. Clearly, the Nernst equation 
cannot fit this case so it must be nonideal. 



Figure 5. Measured and Nernst-predicted open-circuit poten- 
tial. 



Figure 6 . Measured and Redlich-Kister-predicted open- 
circuit potential. 

N p = 12 and N n = 0. Clearly, the modified equation much 
more realistically reflects the tme voltage curve. The iden- 
tified parameters are given in Table 1 and were determined 
using the Nelder-Mead simplex method. The identified pa- 
rameters represent those achieving a least-squares (local) er- 
ror minimum, which is satisfactory for our purposes. 

4.2. Concentration Overpotential 


The activity coefficient terms are related to excess Gibbs free 
energy and can be captured using the Redlich-Kister expan- 
sion (Karthikeyan et al., 2008). We then have 



( 10 ) 


The number of terms iVj in the above expansion and the fitting 
parameters must be identified. Fig. 6 shows the modi- 
fied equation fitted to the measured equilibrium potential for 


When a battery is discharged, the reactions take place at the 
surface of the electrode and this results in a concentration gra- 
dient. The concept is shown in Fig. 7. In the bulk volume, the 
concentration is nearly even, but close to the electrode the 
concentration changes drastically. To accommodate this, one 
may split the total volume into two individual control vol- 
umes (CVs), one for the bulk (with subscript b) and one for 
the surface (with subscript s). 

For the volumes, the concentration of Li ions is computed as 

Qb.i 

Cb,i — 

Vb,i 

Qs,i 

C-s.i — ; 

V s ,i 


(ID 

( 12 ) 
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Table 1. Battery Model Parameters 


Parameter 

Value 

Uo, P 

4.03 V 

A p , o 

-33642.23 J/mol 

A P ,i 

0.11 J/mol 

Ap,2 

23506.89 J/mol 

Ap,3 

-74679.26 J/mol 

Ap,4 

14359.34 J/mol 

Ap, 5 

307849.79 J/mol 

Ap,Q 

85053.13 J/mol 

Ap, 7 

-1075148.06 J/mol 

Ap, 8 

2173.62 J/mol 

Ap,9 

991586.68 J/mol 

Ap, io 

283423.47 J/mol 

Ap, n 

-163020.34 J/mol 

Ap, 12 

-470297.35 J/mol 


0.01 V 

^4n,o 

86.19 J/mol 



Figure 7. Concentration gradient. 


The q variables are described as 


Qs,p 

— ^ app H - Qbs,p 

(17) 

Qb,p 

— Qbs,p ^ app 'Lapp 

(18) 

Qb,n 

— Qbs,n H - 'Lapp ^ app 

(19) 

Qs,n 

^ app Qbs,m 

(20) 


where i app is the applied electric current. Initially, c s ,i = q,j 
and so diffusion is zero. Once discharge begins, ions move 
from the surface layer of the negative electrode to the bulk 
volume of the negative electrode, through the separator to the 
bulk volume of the positive electrode and then on to the sur- 
face layer of the positive electrode. If the applied current is 
larger than the diffusion rate, then a concentration gradient 
will build up between the surface and bulk volumes. Once 
discharge stops, then only diffusion is active, and the concen- 
trations will even out, causing a rise in voltage (because the 
concentration lowers in the surface layer of the positive elec- 
trode, raising the voltage, and the concentration rises in the 
surface layer of the negative electrode, raising the voltage). 

The concentration overpotential is the difference in voltage 
between the two CVs due to the difference in concentration. 
Using the expression for equilibrium potential, we can com- 
pute the potential for the bulk volume and the potential for 
the surface layer; the difference between them is the concen- 
tration overpotential. We can explicitly account for the con- 
centration overpotential simply by using as the expression for 
equilibrium potential, the equilibrium potential of the surface 
layer, i.e., 

Vu,i = Uq-\ — In f - \ + kiNT.jj (21) 

Tlr Y %s,i J 

where x s ,i is computed using 


where, for CV v in electrode i, c v t is the concentration and 
v Vi i is the volume. We define V{ = Vb,i + v s ,i. Note now that 
the following relations hold: 


Qp — Qs, p T Qb,p (13) 

<ln = q S ,n + Qb,n (14) 

Q = Qs, p 4” Qb,p T qs,n T Qb,n- (15) 

As the battery discharges, Li ions must move out of the sur- 
face layer at the negative electrode, through the bulk, and to 
the surface layer at the positive electrode, in order to match 
the flow of electrons. Li ions also move from the bulk CV to 
the surface due to the concentration gradient. The diffusion 
rate from the bulk to the surface is expressed as 

Qbs,i — ~j^j{Cb,i (16) 




Qs,i 
~max ’ 


and 


q 


max 

s,i 


= q 


Vs,i 

Vi 


( 22 ) 


(23) 


The mole fraction in the surface changes faster than the mole 
fraction at the bulk, causing a larger concentration gradient 
for larger applied current and smaller diffusion rates. The 
observed voltage depends only on what is happening in the 
surface layer, not the bulk. When discharge stops, the surface 
layer becomes like the bulk, accounting for the recovery in 
voltage. The battery may then be discharged further since the 
surface layers were replenished with ions (from the bulk) in 
the negative electrode. 


where D is the diffusion constant. 
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4.3. Ohmic Overpotential 

The voltage drops due to the solid-phase ohmic resistances, 
the electrolyte ohmic resistance, and the resistances at the 
current collectors can be treated as constant resistances and 
lumped together: 

V 0 = Vs,p + Vs,n + V e + V CC! p + V cc ,n (24) 

— ‘tappiRs^ ' Rs . n T Re T Rcc,p T Rcc,n ) (25) 

~ Repp Ro ■ (26) 

where V cc is the voltage drop over the current collectors, , 
is the resistance at the current collector for electrode i. 


whereas apparent SOC would be computed based only on the 
surface layer. That is, a battery can be discharged at a given 
rate, and reach the voltage cutoff, i.e., apparent SOC is then 
0. But, once the concentration gradient settles out, the surface 
layer will be partially replenished and the battery can be dis- 
charged further, i.e, apparent SOC increases whereas nominal 
SOC remains the same. 


Nominal (n) and apparent (a) SOC can then be defined using 


SOC n 

SOC a 


Qn 

0.6<7 max 

Qs,n 

0 . 6 ( 7 max »>" ’ 


(31) 

(32) 


4.4. Surface Overpotential 

The overpotentials due to charge transfer resistance and SEI 
kinetics are described by the Butler- Volmer equation. For 
small currents the overpotential is linear with current and for 
larger currents it is proportional to the log of the current, and 
these two characteristics are combined in the Butler- Volmer 
equation. The equation is given as 

) F \ ( aF 

V V p I — exp I — 

(27) 



Ji = Ji o exp 


(1 — a 

rT 


where a is the symmetry factor, Jj is the current density, and 
./,o is the exchange current density. The current densities are 
defined as 

Ji = 4- (28) 

Oi 

Ji 0 = h{ 1 - x Sti ) 0 ‘{x s>i ) 1 ~ a , (29) 

where fcj is a lumped parameter of several constants includ- 
ing a rate coefficient, electrolyte concentration, and maxi- 
mum ion concentration. The voltage drop takes place at the 
electrode surface, so x s/ , is used in computing the exchange 
current density. 

In the case of Li ions, the symmetry factor is 0.5, so the 
Butler- Volmer equation can be simplified and written in terms 
of voltage as 


where g max »." = q m The factor 1/0.6 comes from the 
fact that the mole fraction at the positive electrode cannot go 
below 0.4 (as described in Section 4.1), therefore SOC of 1 
corresponds to the point where q n = 0.6g maXs ". 

4.6. Battery Voltage 

Now that each of the voltage drops in Fig. 4 have been de- 
fined, battery voltage can be expressed as follows. 

V = V UtP - V Uin - K - V vj> - V v>n . (33) 

Voltages in the battery are not observed to change instanta- 
neously, i.e., the voltage changes occur smoothly. When dis- 
charge completes, for example, the voltage rises slowly as 
the surface layers move to the concentrations of the bulk vol- 
umes, as caused by diffusion. In addition to this, there are 
transients associated with V a and the V r/// terms. To take this 
into account in a simple way, we compute voltage using 

v = Vu, p -Vu, n - V' - v^ p - , (34) 

where 

K = (Vo - V')/T 0 (35) 

Kp = C Vv,P - V', p )/Tr,, p (36) 

Kn = (Vn,n - V' tn )/Tr,, ni (37) 

where the r parameters are empirical time constants. 


A. = ^arcsmh (^) . (30) 

4.5. State of Charge 

As discussed in Section 2, the SOC of a battery is convention- 
ally defined to be 1 when the battery is fully charged and 0 
when the battery is fully discharged. In this model, it is anal- 
ogous to the mole fraction x n , but scaled from 0 to 1. There is 
a difference here between nominal SOC and apparent SOC. 
Nominal SOC would be computed based on the combination 
of the bulk and surface layer CVs in the negative electrode. 


The model contains as states x, q StP , q ^ p , qb, n , Qs,n , V' 0 , V p p , 
and V p n . The single model output is V. 

Identified model parameters are given in Table 2. Some pa- 
rameters are defined by the battery dimensions and chemistry. 
The remaining parameters are estimated based on data, as 
with the parameters in Table 1 . 

The measured and predicted discharge curves for a constant 
2 A discharge are shown in Fig. 8. The model fits very well 
to the measured values. The accuracy towards the end of dis- 
charge is most sensitive to the Redlich-Kister parameters, the 
diffusion constant, and the volume of the surface layer. The 
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Table 2. Battery Model Parameters 


Parameter 

Value 

q mdx 

1.32 x 10 4 C 

R 

8.314 J/mol/K 

T 

292 K 

F 

96487 C/mol 

n 

1 

D 

7.0 x 10 6 mol s/C/m 3 

To 

10 s 

a 

0.5 

Ro 

0.085 O 

S P 

2 x 10" 4 nr 

kp 

2 x 10 4 A/m 2 

Vs,p 

2 x 10 _G m 3 

V b,p 

2 x 10 -5 m 3 

T V ,P 

90s 

s n 

2 x 10" 4 m 2 

kn 

2 x 10 4 A/m 2 

Vs,n 

2 x 10 _G m 3 

Vb,n 

2 x 10" 5 m 3 

T~r),n 

90s 



Time (s) 

Figure 8. Measured and predicted 2 A discharge curves. 

predicted end-of-discharge time can be finely tuned most eas- 
ily by adjusting the time at which the amount of Li ions in the 
surface layer of the negative electrode hits zero and the sur- 
face layer of the postive electrode becomes full, because this 
is when the equilibrium voltage of these layers nears a com- 
bined 2.5 V. This is accomplished by adjusting the diffusion 
constant. 

Model validation for a variable loading scenario is shown in 
Fig. 9. As shown in Fig. 9a, the load changes every 2 minutes, 
resulting in corresponding changes in voltage. Fig. 9b shows 
that the voltage predictions are fairly accurate in response to 
changes in load. Some errors are still present that may possi- 
bly be accounted for by including temperature effects. 

4.7. Battery Aging 

As described in Section 2, battery aging manifests in two ma- 
jor ways. The first is a loss of capacity due to parasitic and 




(b) Measured and predicted variable loading discharge curves. 
Figure 9. Model validation for variable loading. 

side reactions that result in a loss of active (mobile) Li ions. 
The second is an increase in internal resistance due to SEI 
layer growth and other factors. Capacity is typically mea- 
sured with respect to a reference current and the associated 
EOD point as defined by a voltage cutoff. A decreased ca- 
pacity will result in the voltage hitting the cutoff earlier. An 
increase in resistance will lower the voltage overall, so also 
result in hitting the cutoff earler. Therefore both these aging 
mechanisms result in a loss of measured capacity. 

Fig. 10 shows 2 A discharge curves with the battery at dif- 
ferent ages. The arrows in the figure illustrate the direction 
the curves move as the battery ages. The EOD point moves 
earlier in time due to diminished capacity. The voltage drops 
down during discharge due to increased resistance. Steady- 
state voltage after discharge increases because the battery is 
not discharged as fully due to the increased resistance and the 
battery reaching EOD before then. 

In the battery model, the total available charge in the battery is 
represented through g max . Therefore, the loss of active mate- 
rial can be represented in the model through a change in g max . 
Fig. 1 la shows how the discharge curve changes as q max is de- 
creased by 1% with each new discharge. We can see a slight 
decrease in voltage, which is due to the Butler- Volmer term. 
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Figure 10. Measured 2 A discharge curves with battery aging. 


which is dependent on mole fraction that is computed using 
q max . The EOD point moves earlier due to the discharge point 
being reached with less charge extracted. Steady-state volt- 
age after discharge increases since the cutoff is reached with 
less relative charge extracted. 

One way to validate this is by looking at the equilbrium po- 
tential plots as the battery ages, because the resistances are 
negligible. Fig. 12 shows the measured equilibrium poten- 
tial of the battery (Vu p — Vjj_ n ) after 10 discharge cycles and 
after 30 discharge cycles. The curves are exactly the same, 
only the time scale changes; there is a 20 minute difference in 
EOD. Since equilibrium potential is a function only of mole 
fraction, and mole fraction is computed as q/q mllx , then de- 
creasing q msx would change the time scale in this way. 

Increase in internal resistance can be captured through an in- 
crease in the R a parameter. Fig. 1 lb shows how the discharge 
curve changes as R 0 is increased by 5% with each new dis- 
charge. The voltage drops, but the change in EOD is not 
significant, because the dynamics near EOD are dominated 
mainly by the equilibrium potential contribution with some 
contribution from the Butler- Volmer dynamics. The change 
in charge extracted is thereby also not significant so changes 
in steady-state voltage after discharge are very small. 

From Figs. 11a and 11b, it is clear that changes in both 
g max and R a are required to capture how the discharge curve 
changes with aging. Fig. 11c shows the combined effects, 
with q max decreasing by 1% and R a increasing by 5% with 
each new discharge. The qualitative changes observed in 
Fig. 10 are now reproduced. 

5. Application to Prognostics 

With an accurate model and known future inputs to a system, 
prognostics should in turn be very accurate. In this section 
we demonstrate battery prognostics with the new model de- 





Figure 11. Simulated 2 A discharge curves with aging. 
Curves for later cycles are shown in lighter colors. 

veloped in this paper. We use the architecture described in 
Section 3. As an estimation algorithm, we use the unscented 
Kalman filter (UKF) with the battery model; see (Julier & 
Uhlmann, 1997, 2004) for details on the filter and (Daigle, 
Saha, & Goebel, 2012; Daigle, Saxena, & Goebel, 2012) 
for its application to prognostics. The UKF operates on a 
set of determinstically selected samples, called sigma points, 
that are used to represent the joint state-parameter distribution 
p(x(k), 9(k)\y{k 0 :k)). 

For the prediction algorithm, we perform a simple simulation 



9 


Annual Conference of the Prognostics and Health Management Society 2013 



Figure 12. Shifting of equilibrium potentials due to aging. 



as described in (Daigle & Goebel, 2013). Each sigma point 
is simulated forward using the model until EOD is reached; 
from the corresponding EODs for each sigma point we can 
construct the EOD distribution. In this work, we assume that 
the future inputs ( i app ) are known, so the only uncertainty 
present in the prediction is that related to the model. We use 
Veod = 2.6 V as the voltage cutoff defining EOD. 

As performance metrics, we use percent root mean square 
error (PRMSE) for estimation accuracy, and for prediction 
we use relative accuracy (RA) and relative standard deviation 
(RSD). We plot results in a- A plots, where a (e.g. 10%) de- 
fines an accuracy cone around the ground truth, and A is a 
time point (Saxena, Celaya, Saha, Saha, & Goebel, 2010). 

Fig. 13 shows prognosis results from a 2 A discharge, assum- 
ing the future inputs ( i app ) are known. As shown in Fig. 13a, 
the UKF tracks very well the measured voltage, since the un- 
derlying model is very accurate. Prediction results are, in 
turn, also very accurate and with very little uncertainty (the 
only uncertainty is related to the state estimation uncertainty), 
as shown in Fig. 13b. Overall, RA averages over all predic- 
tion points (every 100 s) to 99.38%, and RSD to 0.75%. Since 
the spread is so small, the 5%-25% and 75%-95% ranges 
cannot be seen on the plot. 

Fig. 14 shows prognosis results from a 3 A discharge. In the 
open loop, the model slightly underestimates EOD, as shown 
in Fig. 14b. Even though the UKF tracks well, once a pre- 
diction is made the model is in the open loop, and since the 
model predicts an earlier EOD for the given load there is a 
bias in the predictions. The UKF helps keep the state esti- 
mate accurate and so reduces the error compared to using the 
model itself. 

Fig. 15 shows prognosis results from a variable loading dis- 
charge, with EOD being defined in this case as 3.35 V. In the 
open loop, the model slightly overestimates EOD, as shown 
in Fig. 15b. The UKF tracks well and corrects for some of 
the model errors (see Fig. 9b). The model is more accurate 
at predicting EOD when defined towards the end of a dis- 
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Figure 13. Prognostics results for 2 A discharge. 


charge rather than in the middle of a discharge, so if the volt- 
age threshold is defined to be higher, the predictions will be 
less accurate. In this case RA averages to 88.41%, and RSD 
to 2.13%. RSD is larger because the UKF needs to account 
for the additional model uncertainty since the model is less 
accurate in variable-loading scenarios. 

Additional results are shown in Table 3. PRMSE is very low 
as the UKF is able to track the voltage well in all cases. The 
predictions are also very accurate, with average RA above 
90%, and precise, with RSD around 1% or lower. Although 
overall very good, the results are best around the 2 A dis- 
charge, suggesting that further model improvements (e.g., in- 
cluding temperature effects) may be necessary. 

It is also worth mentioning that the model is not computation- 
ally intensive. The model is implemented with a fixed-step 
discrete-time simulation with a step size of 1 s. It takes on 
the order of 1 x 10 -4 s to simulate 1 s of real time. 

6. Conclusions 

In this paper, we presented a new model for battery prognos- 
tics based on the underlying physics. Using equations based 
in electrochemistry, we developed a model that is both accu- 
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Figure 14. Prognostics results for 3 A discharge. 


Figure 15. Prognostics results for variable loading discharge. 


Table 3. Prognostics Performance 


f'app 

PRMSE RA 

RMADbx/l 

1.0 

0.19 

92.77 

1.07 

1.5 

0.17 

96.02 

0.88 

2.0 

0.17 

99.38 

0.75 

2.5 

0.26 

97.75 

0.82 

3.0 

0.41 

96.08 

0.92 


rate and efficient. As a result, prognostics results for EOD 
prediction are very accurate, with the uncertainty associated 
with the model remaining very small. Such a model has uses 
in other areas as well, such as control. 

The model described here can also be applied to prognostics 
of battery packs. There are two approaches here: either each 
battery is modeled individually, or the batteries of the pack 
are lumped into an equivalent single -cell model. In the latter 
case, the prognostics framework as described here applies di- 
rectly. In the former case, Te must be redefined on a system- 
level (Daigle, Bregon, & Roychoudhury, 2012), such that the 
battery pack is said to be at EOD (or EOL) when any single 
battery in the pack is at EOD (or EOL). 

The current model is limited in several ways that provide the 


basis for future work. For one, we have not described how 
the internal battery temperature changes over time. This is 
important because many of the potentials in the battery are 
functions of temperature. Further, diffusion rates should be 
impacted by temperature. Although we have described how 
aging can be represented within the model, we do not yet have 
descriptions of how these age-related parameters change over 
time, i.e., we require damage progression equations. Much 
work has already been done in this area but not with prog- 
nostics requirements and applications in mind. With such 
models, accurate end-of-life predictions can be made. Upon 
adding these missing elements, application of this modeling 
framework to different battery chemistries is also of interest. 
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